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Abstract 

Prudent walks are self-avoiding walks on a lattice which never step into the direc- 
tion of an already occupied vertex. We study the closed version of these walks, called 
prudent polygons, where the last vertex of the walk is adjacent to its first one. More 
precisely, we give the half-perimeter generating functions of two subclasses of pru- 
dent polygons on the square lattice, which turn out to be algebraic and non-D-finite, 
respectively. 

1 Introduction 

The enumeration of self-avoiding walks (SAW) and polygons (SAP) on a lattice by their 
number of steps [TT] is a long standing problem in combinatorics. Extrapolation of series 
data from exact enumeration has led to high precision estimates of the exponential growth 
rate and subexponential corrections but an exact solution of either problem (i.e. finding 
the generating function, see below) seems out of reach. Rechnitzer [16] has shown that 
the anisotropic generating function of SAPs on the square lattice is not D-finite. A (pos- 
sibly multivariate) function /(z) is D-finite, if the vector space over C(z) spanned by its 
derivatives is finite dimensional. In the univariate case this means that / is a solution of a 
homogeneous linear ordinary differential equation with polynomial coefficients. At present 
one tries to find solvable subclasses with large exponential growth rates. This approach 
is particularly successful in two dimensions. We will restrict to the square lattice in our 
paper. 

One promising class are so-called prudent walks (PW) [5l [H], which are SAWs never 
stepping towards an already occupied vertex. Note that a general prudent walk is not 
reversible, i.e. the walk traversed backwards from its terminal vertex to its initial ver- 
tex may not be prudent. Since SAWs are counted modulo translation, we may choose 
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the initial vertex of a PW to be the origin (0,0). The full problem of PW is unsolved, 
but recently Bousquet- Melon [2] succeeded in enumerating a substantial subclass. We 
adopt the terminology of her paper and use the same methods to obtain the generating 
functions for the corresponding polygon models defined below. Every nearest neighbour 
walk on the square lattice has a minimal bounding rectangle containing it, referred to 
as the box of the walk. It is easy to see that each unit step of a prudent walk ends on 
the boundary of its current box. (This is not a characterisation of PWs, e.g. the walk 
(0,0) ^ (0,1) ^ (1,1) ^ (2,1) ^ (2,0) ^ (1,0) is not prudent.) This property allows 
the definition of the following subclasses. Call a PW one-sided, if every step ends on the 
top side and two-sided, if every step ends on top or on the right side of the current box. 
Similarly, a PW is referred to as three-sided if every step ends on the left, top or the right 
side of the current box and additionally each left step and each right step that ends on the 
bottom side of the current box inflates the box. 

Remark, i) As soon as the width of the box of a PW is greater than one, the latter 
additional condition is redundant. It rules out certain configurations which can occur 
only if the box width is equal to one, namely "downward zig-zags" of width one, e.g. 
...^(1,0)^ (1, -1) ^ (0, -1) ^ (0, -2) ^ (1, -2) ^ (1, -3) ^ . . . which needlessly 
complicate the computations below. 

a) Duchi [5j introduced two-sided and three-sided PWs as type-1 and type-2 PWs, respec- 
tively. In [H [7j the authors also employ her notation. 

Explicit expressions for the generating functions of one-, two- and three-sided PWs have 
been found so far confirming the data obtained in [U [7] by computer enumeration. The 
first class consists of partially directed walks and has a rational generating function. The 
second class was shown to have an algebraic generating function by Duchi [5] and recently 
in [2] the third class was solved and the generating function was found not to be D-finite. 

Guttmann [U [71 [9] proposed to study the polygon version of the problem, meaning 
walks, whose last vertex is adjacent to the starting vertex. We exclude single edges from 
this definition. As above, the property of being prudent demands a starting vertex and a 
terminal vertex. So prudent polygons are rooted polygons with a directed root edge. Note 
further that a prudent polygon (PP) which ends, say, to the right of the origin (i.e. in the 
vertex (1, 0)) may never step right of the line x = 1, and furthermore if the walk hits that 
line it has to head directly to the vertex (1,0). So prudent polygons are directed in the 
sense that they contain a corner of their box. Moreover, a fc-sided PP can be interpreted 
as a (fc — l)-sided PW confined in a half-plane, see also Section |6} In this paper we deal 
with the polygon versions of the two-sided and three-sided walks, referred to as two-sided 
and three-sided PPs. Enumeration of one-sided PPs is trivial, since these are simply rows 
of unit cells. We give explicit expressions for the half-perimeter generating functions of 
two-sided PPs and three-sided PPs and show that the latter is not D-finite, as expected 
on numerical grounds [H l^- To our knowledge three-sided PPs are the first exactly solved 
polygon model with a non-D-finite half-perimeter generating function. Enumeration of the 
full class of PPs remains an open problem, as for the walk case. 

Outline: In Section [2] we give functional equations for the generating functions which are 
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based on decompositions of the classes in question, in Section [3] we solve those by the kernel 
method [H 121 [3l [12] and in Section |4] we study the analytic behaviour of the generating 
functions of two-sided and three-sided prudent polygons. Section [5] is dedicated to the 
random generation of PPs. 

2 Functional equations 

In combinatorial enumeration of objects from a class V (say PPs) with respect to counting 
parameters ci, . . . , Cn : P — > Z>o (say perimeter, area etc.) the (multivariate) power series 



P (Xi, . . . , Xn) — ^ ^ X-y^ 



is called the generating function of V. In the following we will count prudent polygons by 
half-perimeter and other, so-called catalytic counting parameters. The variables in the gen- 
erating function marking the latter are called catalytic variables. Their introduction allows 
us to translate certain combinatorial decompositions into non-trivial functional equations 
for the associated generating functions [8]. Furthermore, we identify a PP (a "closed" 
PW) with the collection of unit cells it encloses and build larger PPs from smaller ones by 
attaching unit cells in a prudent fashion, i.e. the new boundary walk with the same initial 
vertex remains prudent. 

A two-sided prudent polygon either ends at the vertex above the origin or at the vertex 
to the right of it. This partitions two-sided PPs into two subsets, which can be transferred 
into each other by the reflection in the diagonal x = ?/. So it suffices to enumerate prudent 
polygons ending on the top of their box. Here two cases occur, namely the "degenerate" 
case when the first steps of the walk are left. The resulting PP is simply a row of unit cells 
pointing to the left. These have a half-perimeter generating function 1^+1^+. . . = t^ / (1— t). 
In the "generic case" such a PP is a bar graph turned upside down, i.e. a column convex 
polyomino containing the top side of its bounding box, cf. Figure [l} Denote by B{t,u,w) 



n_r~L 

Figure 1: Degenerate (left) and generic 2-sided PPs ending on the top of the box 



the generating function of bar graphs counted by half-perimeter, width and height of the 
rightmost column (catalytic parameter), marked by t, u and w respectively. Here w is the 
catalytic variable. The width parameter is not a catalytic parameter. However, it will be 
important in the study of three-sided PPs. We follow the lines of [Ij. 
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Lemma 2.1. The generating function B{t, u, w) of bar graphs satisfies the functional equa- 
tion 



B {t, u,w) = u 



t^w wt{B {t,u,l) - B {t,u,w)) B {t,u,w)t'^w 



1 - wt 



+ 



1 — w 



1 - wt 



(2.1) 



Proof. A bar graph is either a single column, or it is obtained by attaching a new column 
to the right side of a bar graph. The decomposition is sketched in Figure [2] 



j\ + ^ + ^ 



Figure 2: Illustration of the decomposition underlying functional equation (2.1) 



Single columns of height > 1 contribute ut^w/{l — wt) to the generating function. The 
polygons obtained by adding a column which is shorter than or equal to the old rightmost 
column contribute the second summand. This is seen as follows. A polygon of half- 
perimeter n, width k and rightmost column height / contributing t'^u'^w^ to B{t, u, w) gives 
rise to / polygons whose rightmost column is shorter or equal. Their contribution sums up 
to 

, , . t'^w'^l' - t"'U^w'' 

^u^w^ = tuw . (2.2) 



tu ^ t"i 



w 
i=i 

Summing this over all polygons gives the second summand. The third summand corre- 
sponds to adding a larger column. To this end duplicate the rightmost column and attach a 
non-empty column below the so obtained new rightmost column. A so obtained bar graph 
can be viewed as an ordered pair of a bar graph and a column. The generating function 
of those pairs is the third summand of the rhs. This finishes the proof. □ 

The walk constituting the boundary of a three-sided PP has (0, 0) as its initial vertex and 
(1,0) or (—1,0) or (0,1) as its terminal vertex. Those walks with terminal vertex (0,1) 
may not step above the line y = 1 and they have to move directly to the vertex (0, 1) as 
soon as they step upon that line. This leads to two sorts of bar graphs either rooted on 
their left or on their right side, see Figure |3| 



Figure 3: Three-sided PPs with terminal vertex (0, 1) are bar graphs 



So only those three-sided PPs are of further interest, which end in (1,0) or (—1,0). Both 
classes are transformed into each other by a reflection in the line x = 0. We study those 
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ending to the right of the origin in the vertex (1,0). Again a degenerate and a generic 
case are distinguished, according to whether such a PP reaches its terminal vertex from 
below via the vertex (1, —1) ("counterclockwise around the origin") or from above, via the 
vertex (1,1) ("clockwise"). In the degenerate case we simply obtain a single column. In 
the generic case, a (possibly empty) sequence of initial down steps is followed by a left 
step. So denote by R(t, u, w) the generating function of the generic three-sided PPs ending 
in the vertex (1,0) counted by half-perimeter, the length of the top row and the distance 
of the top left corner of the top row and the top left corner of the box, marked by t, u, and 
w, respectively, cf. Figure |4j Here both u and w are catalytic variables. 




Figure 4: Degenerate and generic three-sided PPs, catalytic variables 



Lemma 2.2. The generating function R{t,u,w) of generic three-sided PPs satisfies the 
functional equation 



ut (R [t, w,w) — R (t, u, w)) 



R(t, u, w) = ut (B (t, u)+t) + 

w ~ u 

ut'^ {R {t, u,w) — R (t, u, ut)) 



(2.3) 



+ 



w — ut 



+ R {t, u, ut) ut {B (t, u) + t) 



where B{t, u) := B(t, u,l) is the generating function of bar graphs counted by half-perimeter 
and width. 

Proof. The decomposition we use is sketched in Figure |5} 



i_r 



7^rn + 




+ 



[ri| 



Figure 5: Illustration of the decomposition underlying functional equation (2.3) 



The polygons in question contain the top right corner of their box. This corner is some 
point (1, y). If y = 1, then the PP is either the unit square containing (0, 0) and (1, 1) or a 
bar graph as above with that unit square glued to the right. This yields the first summand. 
A PP with y > 1 is obtained in one of the following three ways from a PP with top right 
corner {l,y — 1). The first is to add a new row on top, which is shorter than or equal to the 
original top row. A similar computation as in (2.2) (with some additional book keeping 
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on w) yields the second summand. The second way to obtain a larger PP from a smaller 
one is by adding a new row on top, which is longer than the original top row, but does not 



inflate the box to the left. Again a treatment similar to the computation in (2.2) yields 



the third summand. The third way to extend a PP is to add a row on top of length equal 
to the width of the box plus one and possibly an arbitrary bar graph. This finally yields 
the fourth summand and the functional equation is complete. □ 

Remark. As in the case of general SAPs [lOj we can define a concatenation of two three- 
sided PPs. Roughly speaking, the one PP can be enlargened by inserting the other one at 
the top corner of the leftmost column, see Figure |6l 




Figure 6: Concatenating two 3-sided PPs 



The numbers pp^^^ of three-sided PPs hence satisfy pp^^'^^^ > pPa™^ ■ PP^^^ ■ This implies 
the existence of a connective constant /3, i.e. a representation pp^^^ = exp(/3m -|- o(m)). 
The precise value for /3 and the subexponential corrections are given in Section |4j The 
converse inequality holds for prudent walks, since breaking an m -|- n step PW after m 
steps leaves one with a pair of prudent walks of respective lengths m and n. 

We now turn to unrestricted PPs. They can be partitioned into eight subclasses accord- 
ing to their end point being (1, 0), (0, 1), (—1, 0) or (0, —1) and their orientation (clockwise 
or counterclockwise around the origin). All eight classes can be transformed into each other 
by symmetry operations of the square. Hence it suffices to enumerate those PPs ending 
in (1,0) which reach their endpoint via the vertex (1, 1) (clockwise). We denote this class 
by J-" and by F{u,w,x) := F(t,u,w,x) its generating function. G{u,w,x) := G(t,u,w,x) 
and H{u, w, x) := H(t, u, w, x) are defined as the generating functions of the two auxiliary 
subclasses Q ^Ti of specified below. We have the following functional equation. 

Lemma 2.3. The power series F , G and H satisfy a system of functional equations. For 
X = F,G, H the single equations are of the form 

, tux (X (w,w,x) — X (u,w,x)) 

X[u,w,x) = 

w~u .2.4) 
t'^ux (X (u,w,x) — X (u,ut,x)) 

H 1 h Ix{u,w,x), 

w — ut 
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where the formal power series Ix{u,w,x) := Ix(t,u,w,x) is equal to 



G{x, X, u) if X=F, 

Ix{u, w, x) := { t^ux + t^uxF{x, xt, w) + xH{x, x, u) if X=G, 
t'^uxw~^G{x,xt^w) ii X=H. 

Proof. The proof relies on a decomposition similar to that of three-sided PPs. The PPs 
of the class J-" all contain the top right corner of their box. In the generating function 
F[u,w,x) of F the variable u marks the length of the top row, w marks the distance of 
the top left corner of the top row to that of the box and x marks the height of the box. 

We define the classes Q and H. Q consists of the unit square together with those PPs 
in F which are obtained by attaching a piece (a collection of unit cells) on top of a given 
PP in F, such that the top side of the box is shifted by one unit, the left side by at least 
one unit and the bottom side by an arbitrary number of units, see Figure [7) 





Figure 7: Illustrations of the classes F, Q and H, catalytic variables 



Note, that the polygons in Q contain the top left corner of their box. The catalytic variable 
u in the generating function G{u, w, x) marks the length of the leftmost column, w marks 
the distance of the lower left corner of the leftmost column to the bottom left corner of 
the box and x marks the width of the box. The class H is the subclass of Q obtained 
by glueing a piece to the left of the leftmost column and thereby shifting the left side of 
the box by exactly one unit and the bottom side by at least one unit, see Figure [7} The 
variable u in H{u, w, x) marks the length of the bottom-most row, w marks the distance 
minus one of the lower right corner of that row to the bottom right corner of the box. The 
variable x marks the height of the box. 

Now the functional equations are derived similarly to the three-sided case. Every PP 
in F can be extended by adding a new row on top which is shorter than or equal to the 
old top row or longer than the old top row, but does not inflate the box to the left. These 
two operations contribute the first and the second summand in the equation for X = F, 
as in the proof of equation (2.3). Inflating the box to the left yields a PP in Q, explaining 
the expression for Ip. 

As for the functional equation G{u, w, x) of the class Q the first two summands on the 
rhs correspond to adding a new column to the left, the expression xH{x,x,u) to adding 
a piece which shifts the bottom boundary of the box, in an analogous fashion as above. 
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The unit square contributes t'^ux, the term t'^uxF{x,xt,w) corresponds to the "minimal" 
polygons in Q obtained by adding a top row on an arbitrary PP of length equal to the 
width of the box plus one. 

The minimal PPs in "H are those obtained by extending a PP in Q adding a column to 
the left of length equal to the height of the box plus one. This explains the term for Ih 
The rest of the rhs corresponds to adding a new bottom row. □ 



3 Solution by the kernel method 

The following result has already been obtained in fT5] as the solution of an algebraic 
equation which arises from a different decomposition of the class. We derive it here for 
completeness and to recall the "classic" kernel method as applied in [1]. 

Theorem 3.1. The generating function B{t,u) := B{t,u,l) of bar graphs counted by 
half-perimeter and width is equal to 



l-t-u{l + t)t- v/t^(l - tyu^ - 2t (1 - t^) u+{l^ 
Bit,u) = . (3.1) 



(3.2) 



Proof. The functional equation |2.1| is equivalent to 

= [t'^uw{l -w)- uwt{l - wt) - (1 -«;)(! - wt)) B{t, u, w) 

+tuw{l - wt)B{t, u, 1) + t'^uw{l - w) 

The kernel equation 

= (t'^uw{l — w) — uwt{l - wt) - (1 - w){l - wt)) 

is a quadratic equation in the catalytic variable w and has the following unique power 
series solution q(t,u) 

, , 1 + (1 - M)t + Ut^ - JtHl - t)2u2 _ 2t (1 - t2) M + (1 - t)2 ^ ^ 

g(t, u) = ^ ^ ^ ^ J ^ ^ -. (3.3) 



Upon substituting w = q{t,u) into Eq. (3.2), the terms with B{t,u,w) are cancelled and 
we can solve for B(t,u, 1), which leads to (3.1). □ 

Remark. In principle, B{t,u,w) can also be computed, by substituting the result for 
B{t,u,l) into Eq. 

By setting u = w = 1 in the bar graph generating function, adding the contribution 
of the degenerate two-sided PPs and multiplication by 2, we obtain the following result so 
far conjectured by series extrapolation from exact enumeration data [1]. 
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Corollary 3.2. The generating function of two-sided prudent polygons is equal to 
1 /I -3t + t2 + 3t3 



v/(l -t){l-3t-t^- t3) 



t \ 1-t 

= A + 6 + 12 + 28 + 72 + 196 + 552 z^ + 1590 
+ 4656 z^^ + 13812 z^^ + 41412 + 125286 z^^ + 381976 



(3.4) 



+ 1172440 z'^ + 3620024 z'" + 11235830 z" + 35036928 z 
+ 109715014 z^^ + 344863872 z^^ + . . . 



,18 



Now we turn to the three-sided case. Note that the sum of the catalytic counting 
parameters, namely the length of the top row and the distance of its top left corner to the 
top left corner of the box, is equal to the width of the polygon. We have the following 
result for the generic three-sided PPs ending on the right. It is derived in a similar way as 
the corresponding result on PWs in [2]. 



Theorem 3.3. The functional equation 2^ has a unique power series solution. For the 
generating function R{t, w, w) of generic three-sided prudent polygons ending on the right 
and counted by half-perimeter and width we have an explicit expression as a an infinite 
sum of formal power series 



k-l 



R(t, w, w) 



j:L{{tqrw)l[K{{tqy 

k>0 j=0 



W 



Here 



q:=qit,l) 



t'^ + 1- Vl -4t + 2t2 + t4 



2t 



with q{t,u) as in (3.3) in the proof of Theorem 3.1. K and L are given by 



and 



Kiw] 



L(w) 



(1 - t) g - 1 - ((1 - t + t^) g - 1) (E (t, qw) + t) w 
l-t{l + t)q-{t{l-t-t'^)q + t^){B (t, qw) + t) w 

{I + t^ - {I - 2t + 2t^ + t^) q) (B (t, qw) + t) w 
l-t{l + t)q- {t{l-t-t3)q + t^){B (t, qw) +t)w' 



(3.5) 



(3.6) 



(3.7) 



(3.^ 



where B{t,u) is the generating function of bar graphs as in (3.1). 
Proof. The functional equation ( 2.3[ ) is equivalent to 

= {ut^{w — m) — ut{w — ut) — {w — u){w — ut)) R(t, u, w) 
+ [ut {B{t, u) + t) {w — u){w — ut) — t^u{w — m)) i?(t, M, ut) 

+ut{w — ut)R(t, w, w) 
+ut{w -u){w- ut) {Bit, u) +t). 



(3.9) 
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We first solve the kernel equation 

{ut^{w — u) — ut{w — ut) — {w — u){w — utf) = 
for u and w. The unique power series solutions are U(t, w) = q(t)w resp. W{t, u) 



q(t)tu, 



with q{t) as in (3.6). We substitute w = W{t,u) in Eq. (3.9) and obtain an expression for 
R{t, u, ut) in terms of R (t, qtu, qtu) , namely 

{q-l)R (t, qtu, qtu) + (1 - g) (1 - tq) {B (t, u) + t) u 



R (t, u, ut) 



1 - gt - (1 - g) (1 - tq) {B (t, u)+t)u 



(3.10) 



Substitute this into Eq. (3.9) and setu = U {t, w). This relates R{t, w, w) and R (t, wtq"^, wtq^) 
as follows: 

R{t, w, w) = K{w) ■ R {t, wtq^, wtq^) + L{w), (3.11) 

with K{w) and L{w) as in (3.7) and (3.8). K{w) and L{w) are a formal power series in t, 
which is seen as follows: B {t, qw) is well-defined as a formal power series in t as [t^~\ B{t, u) 
is a polynomial in u of degree at most — 1. Furthermore by the definition of B we see 
B(t,u) = t^u + 0{t^). The denominator is now easily seen to be 1 + 0(t), so both K{w) 
and L{w) are well-defined as formal power series in t. Inspecting the first few coefficients 
we see {1 — t)q — 1 = O {t^) and 1 — {1 — t + 1^) q = O {t^) , so the numerator of K{w) is 
O {t^) . In a similar way the numerator of L{w) is seen to be w ■ O {t^) . Moreover we have 
tq^ = t + O (t^) . So we can iterate Eq. (3.11) and obtain formula (3.5). □ 



Remark, i) We have the following alternative expressions for K{w) and L{w) : 

((1 - g)(l - qt)qwt {B{t, qw) + t) + t\{q - 1)) {q - 1) 



K{w) 



and 



Liw) 



q{l - qt)^ ((1 - q)qwt {B{t, qw) + t) + t) 

{1 - qt){l - qH){q - l)qtw {B{t,qw)+t) 
g(l - gt)2 ((1 - q)qwt {B{t, qw) +t)+t)' 



(3.12) 



(3.13) 



The expressions (3.7) and (3.8) were obtained by expressing powers of q in terms of q, e.g. 

g2 = (^2 + 1) g _ t) /t^^ 
g3 = (^4 ^2t^ -t + l) q-t^ -t) /^^ 

g4 = (tq [t^ + St" - 2t^ + 3t^ -2t + l) -t + t^ - 2t' - t') /t\ 

ii) In principle one could also compute R(t,u,w). To obtain the generating function of all 
three-sided PPs we sum up the contributions of the degenerate PPs and those ending on 
top, multiply by two and obtain 



t' 



1 -t 



B{t,l) + R{t,l,l) 



The first few terms of the series PP^it) are 

PP3 (t) = 6 + 10 ^3 + 24 + 66 1^ + 198 1^ + 628 + 2068 1^ 



7004 r + 24260 1 
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+ 85596 1^^ + 306692 1^^ + 1113204 1^^ + 4085120 1^^ + 15131436 1^^ + 56495170 1 
+ 212377850 1^' + 803094926 + 3052424080 1^^ + 11653580124 1^^ + . . . . 



16 
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4 Analytic properties of the generating functions 



So far we have considered the generating functions in question as formal power series. A 
crude estimate on the number of SAPs of half-perimeter n is 4^" which is the total number 
of all nearest neighbour walks on the square lattice of length 2n. So the series PP2{t) and 
PP^{t) converge at least in the open disc {\t\ < 1/16} and represent analytic functions 
there. This section deals with the analytic properties of these functions. We first discuss 
the analytic structure of the generating function of two-sided PPs. 



Proposition 4.1. The generating function PP2{t), of. (3.4), is algebraic of degree 2, with 
its dominant singularity a square root singularity at t = p, where p is the unique real root 
of the equation 

l-4t + 2t2 + t4 

— = l-3t-t^ -t^ = 0. 4.1 

1-t ^ ' 



With 9 = ^/26+~6^/^ the exact value for p can be written as 



9^-9-8 

p = — = 0.2955977 .... 

36' 

The number pp^^ of two-sided PPs of half-perimeter m is asymptotically 

pp^^'^ A ■ p-"" ■ m-^/^ (m^cx)), 

where 



-37 + llv^)02 + (_i52 + 8^33)0 + 32 

A = = = 0.8548166 .... 

4 v67rp 

Remark: i) The generating function of two-sided prudent walks is algebraic with its 
dominant singularity a simple pole at cf = 0.403 .... Its coefficients are asymptotically 
equal to k ■ where k, = 2.51 . . . , cf. pj. 

a) The asymptotic number of bar graphs as well as staircase polygons (counted by half- 
perimeter) and Dyck paths (by half-length) is of the form k ■ /x" ■ ra"^/^. Furthermore, the 
area random variables in the fixed-perimeter (fixed-length) ensembles of all three models 
are known to converge weakly to the Airy distribution [6l [T71 [19] . 

The analytic structure of PPsit) is far more complicated due to the analytic structure 



of -R(t, 1,1), which is stated in the main result Theorem 4.4, In what follows we make 
frequent use of the following facts about the series q : 

Lemma 4.2. The series q, (1 — t)q — 1, g^t, t(l -|- t)q and t {1 — t — t^) q + t"^ have non- 
negative integer coefficients. For \t\ < p we have the estimates 



\t 


to 


+ 1 


2\t 





< ' ui , \q t\ < 1, 1(1 - t)g - 1| < p, |1 - t(l + t)q\ > p. (4.2) 
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Equality holds if and only if t = p. Furthermore 

,(p, = ^^ = ±. (4.3) 

The singular behaviour of B {t, q {q^t)^^ and B{t^ qw) plays an important role in the study 
of 

Lemma 4.3. For N > the dominant singularity of B {t, q {q^t)^^ is a^, which is the 
unique solution in the interval [0, p) of the equation 

u{t)-q{qH) =-.—^^-q(qH) =0. 

In particular, a := ctq = = 0.2441312 . . . , where r is the unique real root of the polyno- 
mial t^ +2t'^ +3t —2. The sequence {aiy, N > 0} is monotonically increasing and converges 
to p. Furthermore B{t,qw) is analytic in the polydisc {\t\ < p} x {\w\ < ^/p}■ 

Proof. B{t,u) is singular if and only if 

t2(l - t) V - 2t (1 - t^) u + {l-tY = 0. 

The relevant solution u{t) with u{p) = 1 is 

u{t) = ^. 4.4 

B q (g^t)^ j is singular if q {qH)^ = u{t). This equation has a solution in the interval 

(0,p), as u{t) — )■ 1 and q {q^t)^ -> (p^ + l)/2p = 1/^ > 1, for t p. Here u is strictly 
decreasing and q {qH) strictly increasing. We further see that aN converges to p, as for 
arbitrary fixed t with < t < p we can chose N sufficiently large, such that u{t) > q {q^t)^ , 
see Lemma 4.2[ So aj^ > t, which shows the convergence. Monotonicity follows, as 
q{qH) < q{q^t) for t G (0,p). All these singularities are square root singularities, as 
the expressions under the root are analytic in |t| < p. B(t,qw) is singular, ii w = u(t)/q 
and hence 

I I I^WI ^ r ( \ r 
V"\ = > Vpu{p) = VP, 

with equality if and only if t = p. So there is no singularity inside the polydisc. □ 
Now we are ready to state the main result, which is proven in the subsequent lemmas. 

Theorem 4.4. The function R{t, 1,1) is analytic in the disc {\t\ < a} with its unique 
dominant singularity a square root singularity at a. Moreover it is meromorphic in the slit 
disc 

Da,p = {\t\ <p}\[fT,p), 

and it has infinitely many square root singularities in the set {<Jn, = 0, 1, 2, . . .} . In 
particular, R{t, 1,1) is not D- finite. 
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Remark, i) The number ppg'"^ of three-sided PPs of half perimeter m is asymptotically 
equal to k ■ a~"^ ■ m~^/^ for some positive constant k. In particular, two-sided PPs are 
exponentially rare among three-sided PPs. 

ii) The generating function of three-sided prudent walks has its dominant singularity a 
simple pole at a = 0.403 . . . , as in the two-sided case. It is meromorphic in some larger 
disc of radius p = \/2 — 1 with infinitely many simple poles in the intervall [a, p). Its 
coefficients grow like k, ■ for some k > [2]. 

Possible singularities of R{t, 1,1) in D^ p are zeroes of the denominators of K{w) and 



L{w), places, where the representation (3.5) diverges, and square root singularities of 
B (^t, q (q^t)^^ ■ Now we investigate the analytic properties of the single summands in 
the representation (3.5). 



Lemma 4.5. 1. K (^{q'^t)^^ and L (^{q^t)^^ are analytic in {\t\ < <7j\f}. 

2. K {^q^t)^ vo^ and L (^{q^t)^ are analytic in {\t\ < p} x {\w\ < y/p}- 
Proof. With the above definition of u{t) and a short computation we obtain the estimate 

B{t,q{qHf)\<Bi\tlu{m = y^\. (4.5) 

The denominator of K{w) and L{w) is 

1 - T{t, w) = l-t{l + t)q- {t{l-t-t^)q + t'^){B {t, qw) + t) w. 

T{t, w) is a power series in t and w with non-negative coefficients and T(0, w) = 0. Hence 
we have the estimate 

T (t, {qHf) < T {q{a,fa,f) < T (a., ^) . 
A computation shows that the function 1 — T {t,u{t)/q{t)) 

l'T(^t,'^^=l-t{l+t)q-{t{l-t-t')q + t') (v^ + t)^ 

has no zeroes in [a, p]. This finishes the proof of the first assertion, as K (^{q'^t)^^ and 

L (^{q'^t)^^ do not have poles inside {\t\ < ctn}. Furthermore, the denominator 1 — 

T ^t, (qH)^ is analytic in the polydisc {\t\ < p} x {\w\ < i/p}, with the only singular 
point (t, w) = (p, y/p) on its boundary. As above we see 

u{py 



<T(|t|,|^|)<T p 



lip) 



T (p, Vp) 



and hence the denominator is non-zero in the domain in question and K (j.q'^t)^ n?j and 
L (^{q'^t)^ are both analytic in the polydisc. □ 
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Lemma 4.6. 1. The series representation (3.5) of R{t, 1,1) is a series of algebraic 
functions, which converges compactly in the slit disc D^p = {\t\ < p} \ [o"? p) to a 
meromorphic function. 

2. Furthermore the corresponding representation of R{t,w,w) converges compactly in 
the polydisc {\t\ < p} x < y/p} to an analytic function. 



3. The Taylor expansion of R{t, w, w) about {t, w) 
p} X {\w\ < y/p}. 



(0,0) converges absolutely in {\t\ < 



Proof. For the first assertion choose < r < p. We look at the disc {\t\ < r}. The term 
independent of w in the numerator of K{w) is strictly less than p for \t\ < r and the 
corresponding term in the denominator is strictly larger than p, see Lemma 4^ So we can 

choose N large such that (Tn > t and K [{qH)^^ < 1 for |t| < r. Split the series at A^. 

The summands for = 0, . . . , — 1 sum up to a function which is meromorphic in the slit 
disc {\t\ < r} \ [a, r]. In the rest of the series take out the common factors to obtain 



N-l 



k~l 



n ((tgT) ((^^^)"") ( WO • 



(4.6) 



j=o 



fc>0 



j=0 



The first product is a meromorphic function in the slit disc. L (^(tq^)^^''^ is easily seen 

to converge uniformly to in \t\ < r as — !■ oo. In \t\ < r all summands are holomorphic 
(see the above discussion) and the sum can be estimated by a geometric series and hence 
converges uniformly in the compact disc {\t\ < r}. By Montel's theorem the limit of the 
sum is again analytic. This finishes the proof for the first assertion. The second assertion 
is proven along the similar lines. By the multivariate version of Montel's theorem [18] 
the limit function is also analytic in the domain in question and thus the third assertion 
follows. □ 

Lemma 4.7. R{t, 1, 1) is singular at infinitely many of the (Tn- Furthermore, R{t, 1, 1) is 
singular at a. 

Proof. Terms singular at (Tat only show up in the summands for k > N. The sum of these 
(4.6 ) is equal to 



Af-l 

n ^ {i'^^y 

j=0 



L({tqY)+K({tqY)R(t, {tqY^\{tq') 



2\N 



(4.7) 



In order to show that the singularity does not cancel, only the term in square brackets 
is of interest. Singular terms show up in the numerators and the common denominator 
of K [{tqY) and L ({tq^) ■ We now manipulate the expressions ( |3.12[ ) and ( |3.13[ ) for 
K{w) and L{w) in order to get rid of singular terms in the denominator, where the factor 



q)qwt {B{t, qw) + t) + t 
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leads to a singularity at ctat for w = {tq^)^ ■ Write 



where 



qwt {Bit, qw) +t) = A{w) - (j){w), 
1 



A{w) 



[1 + t - qw{l + t)t) 



4>{w) = ^7^2(1 _ ty^qwy - 2t (1 - t2) qw + {l- tf. 

Then A {^tq^)^^ is analytic in {\t\ < p}. After multiplication of the numerator and de- 
nominator with (1 — q)A{w) + t + (1 — q)(j){w) there is no more occurence of in the 
denominator. We now have to collect the terms involving (f){w) in the numerators of K{w) 
and L[w). In the numerator of K{w) the terms involving (j){w) sum up to 

P^H0H := t{q - 1)2(1 - qH)(l){w). 

The terms involving (f){w) in the numerator of L{w) sum up to 

Pz.H0H := (1 - qt){l - qH){l - q)t<P{w). 

So the singularity at ctjv can only cancel if 

Pl {{crMq{aNff^ 

In order to prove that this equation can hold for at most finitely many of the ctat, we 
show that for sufficiently close to p the Ihs of eq. (4.8) is strictly decreasing while 
the rhs is strictly increasing. Since {(Tn) is monotonically increasing and converges to p 
this will finish the proof. We first prove the assertion on the rhs. The Taylor expansion 
of R{t, w, w) about (0, 0) has non-negative coefficients and represents R(t, w, w) in the 
polydisc {\t\ < p} X {|w| < y/p} by Lemma 4.6, By the definition of aN and u{t) we have 

((TArg (O-Tv)^)^^^ = U ((Tiv) q (ciAr) (Ttv 

The rhs of the last equation is strictly increasing for sufficiently large N and converges to 
y/p as — )■ oo. The sequence is also strictly increasing by Lemma 4.3 So for large 

is strictly increasing. 



R (o-AT, (crArg(crAr)^)^'^^ , {<JNq{<yNf)^j 



(4.8) 



enough the sequence R ( ctat, [a^q (ctat) j 



[o-Nqio-N) ) 



Now we turn to the Ihs of eq. (4.8|). A computation yields 

Pl 



(cr^qMY) _ 1 _ ar^qicTj,) 



Pk ( {o-NqicTN 



^2^N 



q{(^N) 



which easily seen to be ultimately strictly decreasing. This finishes the proof of Lemma 
14^1 □ 



The Lemmas 4.5, 4.6 and 4.7 together constitute a proof of Theorem 4.4 
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5 Random three-sided prudent polygons 



In P] prudent walks of a given fixed length are generated uniformly at random with a refined 
version of a method proposed in [13]. We briefly describe a version of the method tailored 
to our particular needs. The main ingredient in the sampling procedure are generating 
trees. These are rooted trees with their nodes labelled in such a way that if two nodes 
bear the same label, then the multisets of the labels of their children are the same. In 
this section we present generating trees and sampling procedures for the various classes of 
prudent polygons. 



The decompositions underlying the functional equations (2.1), (2.3) and (2.4) (see also 
Figures |2| [5] and [T]) yield rules according to which a larger PP from the respective class can 
be constructed starting from a smaller one. We refine these building steps such that each 
step increases the half-perimeter by one. The result is a step-by-step procedure which allows 
to generate any PP of half-perimeter m in a unique way, starting from the unit square, 
such that after the kth step we have a PP of half-perimeter k + 2, k = 0,1, . . . m — 2. To put 
it differently, we can organise the polygons in a rooted tree T, with the unit square as the 
root and the polygons of half-perimeter m as the nodes on level m — 2. So a random PP of 
half-perimeter m corresponds to a random path of length m — 2 in that tree starting from 
the root. A polygon a is a child of a polygon vr, if it is obtained by one of the following 
six construction steps. 

1. Attaching a new top row which is shorter than or equal to the current top row, 

2. attaching a unit square to the left side of the current top row, 

3. attaching a new leftmost column which is shorter than or equal to the current leftmost 
column, 

4. attaching a unit square to the bottom side of the leftmost column, 

5. attaching a new bottom-most row which is shorter than or equal to the current 
bottom-most row, and 

6. attaching a unit square to the right side of the bottom- most row. 

Any of these steps, if applicable, increases the half-perimeter by one, see Figure [8] below. 



1. 



Figure 8: The types of steps used to obtain generating trees 



Remark, i) Steps of types 2,4 and 6 are only admissible if the current top row (leftmost 
column, bottom row) is longer than or equal to the second row from top (second column 
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from the left, second row from the bottom). Additionally, a step of type 6 is forbidden, if 
the bottom row is only one unit shorter than the width of the box. 



ii) In the proof of the functional equation (2.3) the steps of types 3 and 4 are encapsulated 
in the "attaching a bar graph to the left" operation. Hence any generic three-sided PP can 
be generated starting from the unit square by using only steps of the first four types. 
in) Building bar graphs only requires steps of type 1 and 2. Here we reflected the bar 
graphs discussed earlier in the line x = y. 

In order to compute the appropriate probabilities according to which each step in the 
random path in the tree is chosen, we associate to each PP a label encoding the admissible 
steps which can be applied to enlarge it. We give a labelling for unrestricted PPs, since 
the labellings of the other classes are obtained as specialisations thereof. The labels are 
five-tuples {a,e,k,l,p) G {B,L,T} x {y,n} x Z>q. The letter a encodes the last building 
step. It is equal to T (top), if the last step was of type 1 (which inflated the box to the 
top), or of type 2 but without inflating the box to the left, a is equal to L (left) if the 
last step was of type 2 and thereby inflating the box to the left, of type 3 or of type 4 but 
without inflating the box. Finally, a is equal to B (bottom) if the last building step was of 
type 4 with an inflation of the box, or of types 5 or 6. 

If a = T, the parameter e G {y,n} (yes/no) indicates if the current top row is longer 
than or equal to the second row from the top, and hence if a step of type 2 is applicable. 
Similarly, if a = L, e decides if a step of type 4 can be performed, i.e. if the leftmost 
column is shorter than or equal to the second but leftmost one. Finally, ii a = B, e decides 
whether a step of type 6 can be performed. 

The parameter k always denotes the length of the top row, and I is either the distance 
of the left end of the top row to the left side of the box or the length of the leftmost 
column or the length of the bottom row, depending on whether a = T or a = L or a = B 
respectively. The unit square receives the label (L,n, 1, 1,0). 



Labels for the generating tree of PPs 


a 


e 


k 


/ 


P 


T 


top row extendable? 


top row 


dist. of box to top row 


height 


L 


left col. extendable? 


top row 


length of leftmost col. 


dist. of box to left col. 


B 


bottom row extendable? 


top row 


length of bottom row 


height 



Remark 

polygons with a 



In the proof of equation (2.4) we introduced two subclasses Q and T-i. The 
L (resp. a 



B) are precisely those in Q (resp. "H). 

The construction steps yield the following rewriting rules for the labels associated with 
general PPs. 



(T,n, fc,/,p) 



(T, n,i,l + k — i,p 
{T,y,kJ,p+l) 



1), 



k-1 



(5.1) 
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{T,j,k,l,p) ^ <^ 



{L,n,k,l,p) < 



{L,j,k,l,p) 



{B,n,k,l,p) < 



T,n,i,l + k - i,p+ 1), 
T,Y,k,l,p + l) 
T,Y,k + lJ- 
L,n,fc + 1), 

T,n,i,k - i,l + p+ 1), 
T,j,k,OJ + p + l) 
L,n,k + + p - i), 
L,Y,k + l,l,p) 

T,n,i,k - i,l + p + 1), 
T,Y,k,0,l+p+l) 
L,n,k + + p - i), 
L,Y,k + 1,1, p) 
L,Y,k,l + l,p - 1) 
B,n,k,l,l + 1) 

T,n,i,k -i,p+ 1), 
T,Y,k,0,p+l) 
L,n,k + l,i,p-i), 
L,Y,k + l,l,0) 
B,n,k,i,p+ 1), 
B,Y,k,l,p+l) 
B, n, k,l,p + 1) 



if / > 1 
if Z = 



J- 1 



if p > 1 
if p = 



{B,Y,k,l,p) < 



k-1 



l-l 



if k-1 > I 
if k -1 = 1 



k-1 



,P 



,1 



1 



if k-1 = I 
ifk-l>l 



(5.2) 



(5.3) 



(5.4) 



(T, n, i,k - i,p + 1), 
{T,Y,k,0,p + l) 
{L,n, k + l,i,p- i), 
{L,Y,k + 1,1,0) 
[B,n,k,i,p+1), 
[B,Y,k,l,p + l) 
{B,n,k,l + l,p) 
{B,Y,k,l + l,p) 

The labelled rooted tree generated according to these rewriting rules with its root labelled 
(L, n, 1, 1, 0) is a generating tree for unrestricted prudent polygons (more precisely, for the 
class J-', cf. Section [2]). This tree is of course isomorphic to the tree T defined above having 
PPs as nodes, simply by replacing each PP by its label. 

As mentioned above, choosing a PP of half-perimeter m uniformly at random is equiv- 
alent to choosing an m — 2-step path starting from the root uniformly at random. This is 
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achieved by picking each step in the path according to an appropriate probabihty which 
in turn can be expressed in terms of extension numbers. If vr is a polygon of half-perimeter 
m — s (a path of length m — s — 2), then EX{7r,s) denotes the number of polygons of 
half-perimeter m which can be reached from vr in s construction steps, or equivalently of 
extensions of length s of the path. Denote by Ch{7r) the set of polygons obtained from 
TT in one step, i.e. the children of vr in the generating tree. Now the right probability to 
choose a E Ch^ir) in the course of our random sampling procedure is equal to 

, , EX(a,s-l) 

^ ' ^ EX{7T,S) 

The numbers EX{7r, s) can be computed recursively, namely 

. X fl ifs = 0, 

EX{7i,s) = \^ 

[22aech{n) EX{a,s-l) otherwise. 

The crucial observation is that EX{n, ■) only depends on the label of vr, which allows an 
efficient computation. 

For unrestricted PPs, in the first m — 2 levels of the tree O(m^) different labels occur 
since none of the parameters exceeds m. It hence takes O(m^) operations to compute the 
all required extension numbers. We have implemented the procedure and computed these 



numbers up to m = 80. See Figure 12 at the end of the paper for some samples. 

Modifications for two-sided PPs. As remarked above, generating two-sided PPs only 
requires steps of types 1 and 2. The only required information for the building procedure 
is the length of the top row and if the top row is extendable. We hence only need labels 
(T, n, k) and (T, y, k) obtained from the labels (T, -, /c, above by leaving the parameters 



I and p unconsidered. The rewriting rule (5.1) can be adapted unchanged (up to deleting 



the last two coordinates) and in (5.2) simply omit the last line (and the "if / > l"-clause 
in the second but last line). The unit square receives the label (T, y, 1). There are 0{m) 
different labels in the first m — 2 levels of the generating tree, and hence 0{im}) extension 



numbers have to be computed. See Figure 10 for some samples of half-perimeter 250. 



Modifications for three-sided PPs. For the generation of three-sided PPs steps 1, 2, 
3 and 4 suffice. For an appropriate labelling we can hence dump down the {B,-,k,l,p) 
labels and use labels (T, n. A;,/), {T,j,k,l), (L,n, /c,Z) and {L,j,k,l) obtained from the 
labels {T,-,k,l,p) and {L,-,k,l,p) by simply discarding the parameter p. The rewriting 



rules (5.1), (5.2) and (5.3) are adapted without change. In the rule (5.4) drop the last line. 
The unit square is labelled with (L,n, 1,1). We have 0{w?) different labels on the first 
m — 2 levels of the tree and hence O(m^) extension numbers have to be computed. See 



Figure 11 for some samples of half-perimeter 250. 



6 Conclusion 

We have solved the class of two-sided and three-sided prudent polygons, the generating 
function being algebraic in the former and non-D-finite in the latter case. The analysis 
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shows that two-sided PPs are exponentially rare among three-sided PPs which is different 
from the corresponding walk models where the growth rates are equal. 
It would be nice to solve the class of general prudent polygons. We expect that the 
involved functional equations require three or more catalytic variables, which is the case 
for the equation found for the walk model. 

Since the exponential growth rates of SAWs and SAPs are known to be equal [10] it is 
also interesting to compare the exponential growth rates of fc-sided PWs and PPs. To that 
end it suffices to study PPs ending in (1,0). As already mentioned in the introduction, a 
fc-sided PP ending in (1,0) may never step right of the line x = 1 and it heads towards the 
vertex (1, 0) as soon as it hits that line for the first time in a point (1, yo). Up to that step 
the boundary walk of that A;-sided PP is genuinely k — 1-sided. This yields an injective 
map sending a fc-sided PP to a A; — 1-sided PW simply by reflecting the segment joining 
(l,yo) and (1,0) in the line y = y^, see Figure [oj We denote the so-obtained subclass of 




Figure 9: Embedding of /c-sided PPs into k — 1-sided PWs 



k — 1-sided PWs by "embedded fc-sided PPs". If we count PPs by full perimeter, their 
exponential growth rates become 1/^^ = 1-83 ... for two-sided PPs and \l \fo = 2.02 . . . 
for three-sided PPs. It is known that the exponential growth rate of PWs is equal to 
1 -|- ^/2 = 2.41 ... in the one-sided case and equal to 2.48 ... in the two- and three-sided 
cases |2J. The latter rate is also expected for unrestricted PWs [U [7j. Consequently, for 
k = 2,3, our results imply that fc-sided PPs are exponentially rare among /c-sided PWs 
and, via embedding, among k — 1-sided PWs. Furthermore, the rate of three-sided PPs is 
even smaller than that of one-sided PWs. This is not surprising looking at the pictures in 
Figure 11, as such a PP roughly consists of two "almost" one-sided PWs, one heading to 
the far left followed by one "almost directed" walk up and to the right (and the closing 
tail). We expect that this heuristic argument also applies in the general case, which is also 
supported by an estimated value of approximately 2.1 < 1 + a/2 for the growth rate of 
general PPs [11[7|. 
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Figure 11: Random 3-sided PPs of half-perimeter 250 
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Figure 12: Random unrestricted PPs of half-perimeter 80 
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